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ABSTRACT 

We have studied the relationship among nuclear radio and X-ray power, Bondi rate and the 
kinetic luminosity of sub-Eddington active galactic nuclear (AGN) jets, as estimated from 
the pdV work done to inflate the cavities and bubbles observed in the hot X-ray emitting 
atmospheres of their host galaxies and clusters. Besides the recently discovered correlation 
between jet kinetic and Bondi power, we show that a clear correlation exists also between 
Eddington-scaled kinetic power and bolometric luminosity, given by: \og{L]^ia/ L-^dd) = 
(0.49±0.07) log(Lboi/iEdd) — (0.78±0.36). The measured slope suggests that these objects 
are in a radiatively inefficient accretion mode, and has been used to put stringent constraints 
on the properties of the accretion flow. Interestingly, we found no statistically significant cor- 
relations between Bondi power and bolometric AGN luminosity, apart from that induced by 
their common dependence on Lkin, thus confirming the idea that most of the accretion power 
emerges from these systems in kinetic form. We have then analyzed the relation between ki- 
netic power and radio core luminosity. Combining the measures of jet power with estimators 
of the un-beamed radio flux of the jet cores as, for example, the so-called 'fundamental plane' 
of active black holes, we are able to determine, in a statistical sense, both the probability dis- 
tribution of the mean jets Lorentz factor, that peaks at Fm ^ 7, and the intrinsic relationship 
between kinetic and radio core luminosity (and thus the jet radiative efficiency), that we es- 
timate as: logLkin = (0.81 ± 0.11) logLR + 11.9l4'4, good agreement with theoretical 
predictions of synchrotron jet models. With the aid of these findings, quantitative assessments 
of kinetic feedback from supermassive black holes in the radio mode (i.e. at low dimension- 
less accretion rates) will be possible based on accurate determinations of the central engine 
properties alone, such as mass, radio core and/or X-ray luminosity. As an example, we sug- 
gest that Sgr A* may follow the same correlations of radio mode AGN, based on its observed 
radiative output as well as on estimates of the accretion rate both at the Bondi radius and in the 
inner flow. If this is the case, the supermassive black hole in the Galactic center is the source 
of ~ 5 X 10"^* ergs s^^ of mechanical power, equivalent to about 1.5 supernovae every 10^ 
years. 

Key words: accretion, accretion disks - black hole physics - galaxies: active - galaxies: 
evolution - quasars: general 



1 INTRODUCTION 

The growth of supermassive black holes (SMBH) through mass 
accretion is accompanied by the release of enormous amounts 
of energy which, if it is not advected directly into the hole 
(see e.g. Narayan & Yi 1995), can be either radiated away, 
as in Quasars and bright Seyferts, or disposed of in kinetic 
form through powerful, collimated outflows or jets, as ob- 
served, for example, in Radio Galaxies. The feedback exerted 
by such a powerful energy release on the surrounding gas 
and stars is imprinted into the observed correlations between 
black hole mass and galaxy properties JGebhardt et al. 20001 
IFerrarese & Merritt 20001 IMarconi & Hunt 2003t . as well as into 
the disturbed morphology of the hot, X-ray emitting atmospheres 



of groups and clusters harboring active SMBH in their centers 
(|B6hringer et al. 1993[ fPabian et al. 20001 IChurazov et al. 20021 
IBirzan et al. 20041 IFabian et al. 2006b . Radiative and kinetic feed- 
back differ not only on physical grounds, in terms of coupling ef- 
ficiency with the ambient gas, but also on evolutionary grounds. 
The luminosity dependent density evolution parameterization of 
the AGN luminosity functions (|Hasinger, Miyaji & Schmidt 2005[ 
[Hopkins, Richards & Hernquist 2007^ implies a delay between the 
epoch of Quasar dominance and that of more sedate Radio Galax- 
ies ( Merloni 2004 ; Merloni & Heinz 2006 1, such that kinetic en- 
ergy feedback plays an increasingly important role in the epoch 
of cluster formation and virialization. 

From the theoretical point of view, the idea that AGN 
were responsible for the additional heating needed to ex- 
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plain the cooling flow riddle (see e.g. Fabian et al. 2001 
and reference therein) in clusters of galaxies has been put 
forward by many in recent years 
Begelman 200T] IChurazov et al. 20011 



to"'"'" 

finey & Tabor 1995 ; 
lurazov et al. 2002 . 

[Sijacki & Springel 2006[ 

Cattaneo & Teyssier 2007) . However, only very recently was 



Dalla Vecchia et al. 20041 



the specific importance of mechanical (as opposed to radiative) 
feedback for the heating of baryons within the deepest dark 
matter potential wells fully acknowledged by semi-analytic 
modelers of cosmological structure formation dCroton et al. 2006] 
IBower et al. 2006b . Within these schemes, because the bright 
quasar population peaks at too early times, a so-called "radio 
mode" of SMBH growth is invoked in order to regulate both cool- 
ing flows in galaxy clusters and the observed sizes and colors of the 
most massive galaxies ( [Springel et al. 2005|lCroton et al. 2006l l. 

Such indirect evidence for a radio mode of AGN activity 
resonates with a recent body of work aimed at understanding 
the physics of low-luminosity AGN. Detailed multi-wavelength 
studies of nearby galaxies have revealed a clear tendency for lower 
luminosity AGN to be more radio loud as the Eddington-scaled 
accretion rate decreases (see Ho 2002, and references therein; 
Nagar, Falcke and Wilson 2005). Scaling relations with black 
hole X-ray binaries (BHXRB) ( [Merloni, Heinz & Pi Matteo 2003 1 



[Falcke, Kording & Markoff 2004[ IChurazov et al. 2005t have 
also helped to identify AGN analogues of low/hard state 
sources, in which the radiative efficiency of the accretion flow 
is low ([Pellegrini 2005a{ [Chiaberge,Capetti&Macchetto 2005 [ 



[Hardcast e, Evans & Croston 2007) and the radio emitting jet 
carries a substantial fraction of the overall power in kinetic form 
( [Kording, Fender & MigUari 2006) . 

Nonetheless, despite such a widespread consensus on the im- 
portance of kinetic energy feedback from AGN for galaxy forma- 
tion, the kinetic luminosity of an AGN remains very difficult to 
estimate reliably. Recent incarnations of structure formation mod- 
els (either numerical or semi-analytic) hinge on the unknown effi- 
ciency with which growing black holes convert accreted rest mass 
into jet power. Constraints on this quantity, such as those recently 
presented in Heinz, Merloni & Schwab (2007) are clearly vital for 
the robustness of the models. 

Observationally, progress has been possible recently, thanks 
to deep X-ray exposures of nearby elliptical galaxies and clusters 
which have allowed the first direct estimates of Lkin by study- 
ing the cavities, bubbles and weak shocks generated by the radio 
emitting jets in the intra-cluster medium (ICM) dBirzan et al. 20041 
lAUen et al. 20061 [Rafferty et al. 2006) . Two main results emerged 
from these studies: first, it appears that AGN are energetically able 
to balance radiative losses from the ICM in the majority of cases 
( [Rafferty et al. 2006); second, at least in the case of a few nearby 
elliptical galaxies, there is an almost linear correlation between the 
jet power and the Bondi accretion rates calculated from the ob- 
served gas temperatures and density profiles in the nuclear regions 
dAllen et al. 2006] |. The normalization of this relation is such that 
a significant fraction of the energy associated with matter entering 
the Bondi radius must be released in the jets. 

Here we address a different, complementary issue. By col- 
lating all available information on the nuclear (AGN) properties 
of the sources for which the kinetic luminosity was estimated, we 
construct a sample of sub-Eddington accreting supermassive black 
holes with unprecedented level of information on both the inner 
(black hole mass) and the outer (mass supply through the Bondi ra- 
dius) boundary condition for the accretion flow, as well as a reliable 
inventory of the energy output, either in radiative or kinetic form. 



We demonstrate here that with the aid of this new set of informa- 
tion, strong constraints can be placed on the accretion properties 
of these objects. In so doing, we derive a new, robust estimator of 
the kinetic power of a SMBH based on its nuclear properties alone, 
namely its mass and instantaneous X-ray and radio (core) luminos- 
ity. 

The structure of the paper is as follows: In section|2|we study 
the relationships among nuclear radio and X-ray power, Bondi rate 
and the kinetic luminosity of the AGN in the sample and show that 
a clear relationship exists between Eddington scaled kinetic power 
and hard X-ray luminosity. We then analyze the relation between 
kinetic power and radio core luminosity (^. In section!?] we dis- 
cuss our results and present a simple coupled accretion-outflow disc 
model which is capable to explain the main features of the observed 
sample. Finally, we summarize our conclusions in ^ 



2 THE RELATIONSHIP BETWEEN KINETIC POWER 
AND NUCLEAR BOLOMETRIC LUMINOSITY 

We have collected from the literature data on the nuclear properties 
of AGN with published measurements of the jet kinetic power, as 
estimated from the pdV work done to inflate the cavities and bub- 
bles observed in the hot X-ray emitting atmospheres of their host 
galaxies and clusters. We begin by considering the samples of Allen 
et al. (2006) (9 sources) and Rafferty et al. (2006) (33 sources). By 
taking into account common objects, there are 38 AGN with kinetic 
power estimated in such a way. Out of those, we consider only the 
ones for which black hole mass could be estimated (either through 
the M ~ a relation or via direct dynamical measurements), that 
amount to 21 objects. Finally, we search in the literature for avail- 
able measures of the nuclear luminosity in the radio (at 5 GHz) 
and in the 2-10 keV band. Only 6/21 do not have such informa- 
tion, and we further notice that 4 out of these 6 are at a distance 
equal to, or larger than, that of the most distant object in our final 
sample (with the exception of Cygnus A). Therefore, there are only 
two objects (NGC 708 and ESO 349-010, both from the Rafferty 
et al. 2006 sample) which have a measure of Lkin, an estimate of 
the black hole mass and are close enough (z < 0.055) to qual- 
ify as sample members, but which have been dropped due to the 
lack of X-ray and radio nuclear luminosities. Such a high level of 
"completeness" of the sample should guarantee us against substan- 
tial selection bias due, for example, to beaming effects in the radio 
band (see sectionjSjbelowfl 

Our selection provided us with a sample of 15 objects, 2 of 
which have only upper limits to their 2-10 keV nuclear X-ray lumi- 
nosity. In four cases, estimates of the jet power are also available 
based on detailed modelling of either radio emission in the jets and 
lobes (Cyg A, Carilli & Barthel 1996; M87, Bicknell & Begelman 
1996; Per A, Fabian et al. 2002), or shock in the IGM induced by 
the expanding cocoons (NGC 4636, Jones et al. 2007). In another 
handful of objects (M84, M87, NGC 4696, NGC 6166) both Allen 
et al. (2006) and Rafferty et al. (2006) report (independent) mea- 
sures of the jet kinetic power, which, incidentally, differ on average 
by almost one order of magnitude. In all these case when more than 



^ It is of course important to notice tliat the overall selection criteria of 
the sample discussed here are indeed heterogeneous. In particular, the role 
played by the requirement of having a bright, hot X-ray emitting atmosphere 
against which detect bubbles and cavities is very hard to assess, as surface 
brightness effects and different exposure times of the original observations 
should all contribute to the measurability of kinetic power 
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Table 1. Main properties of the sample studied 



Object 


D 


Log Lr Log Lx 


Log Mbh.ct 


Log Lkin 


ffeondi 


References 


(1) 


(2) 


(3) (4) 


(5) 


(6) 


(V) 


(8) 



Cyg A" 


247 


41 


.43 


44.22 


9.40'' 


45.41+0;!/ 


45.24^= 


1,2,3,4,5 






JO. 


.oU 




8 on 


AA (11 +0.16 

44.01_(j 26 


/1 1 _i_n HQ 


A 7 8 
0, /,o 


NGC 1275 (Per A)" 


77.1 


40, 


.74 


43.40 


8.64 


44.33lOiI 


44.3 1^= 


6,9,4,10 


NGC 4374 (M84) 


17 


38, 


.43 


40.34 


8.80 


42.59^^5 


43 69+°-^° 


6,11,8,4 


NGC 4472 


17 


36, 


.69 


38.46 


8.90 


42 91+0.14 


43 79+0-25 


6,7,8 


NGC 4486 (M87)" 


17 


38, 


.88 


40.55 


9.48'' 


43.44l0;5 


44 15+0-28 


6,12,13,4,8,14 


NGC 4552 (M89) 


17 


38, 


.23 


39.33 


8.57 


42 20+0-14 

tZ.ZU_Q 21 


43 37+0.22 
tJ.J/_Q 21 


6,7,8 


NGC 4636" 


17 


36, 


.40 


<38.40 


8.20 


+0.11 

tZ.OJ_Q j^g 


42 29+°-2'' 
tz.zy_Q 24 


6,15,8,16 


NGC 4696 


44.9 


39, 


.10 


40.26 


8.60 


42 SQ+0-22 
tz.5y_g 22 


43.40tO-je 


17,4,8 


NGC 5846 


24.6 


36, 


.50 


38.37 


8.59 


41.86«-i« 


42.85«;g 


6,7,8 


NGC 6166 


135.6 


39, 


.95 


40.56 


8.92 


43.82 j:^;^ 


43 49+0-34 


18,19,4,8 


IC 4374 


94.5 


40, 


.27 


41.37 


8.57 


43 30+°'^'' 
tJ.JU_Q 25 


44.37'= 


20,4 


UGC 9799 


152 


40, 


.55 


41.89 


8.58 


44 1 8+0-36 


43.92= 


21,22,4 


3C 218 (Hydra A) 


242 


40, 


.91 


42.17 


8.96 


44 63+0-16 


44.90= 


23,24,4 


3C388 


416 


40, 


.69 


41.69 


9.18 


44 30+0.38 


44.80= 


25,4 



Notes: " Objects for which a measure of the kinetic jet power from modelling of either jets and radio lobe emission or shocks exists; dynamical mass 
measurements; = Bondi power extrapolated from measures of gas temperature and density outside the Bondi radius, assuming ar~^ density profile. 
Col. (1): Name of the object. Col. (2): Adopted distance in Mpc. Col. (3): Logarithm of nuclear (core) luminosity at 5 GHz. Col. (4): Logarithm of the 
intrinsic rest-frame luminosity in the 2-10 keV band. Col. (5) Logarithm of the black hole mass as derived from M — a relation. Col. (6) Logarithm of the 

Kinetic Luminosity; when multiple estimates available, the logarithmic mean is used. Col. (7) Logarithm of the Bondi Power as defined in eq. {T) 
REFERENCES: (1) Sambruna et al. (1999); (2) Young et al. (2002); (3) Tadhunter et al. (2003); (4) Rafferty et al. (2006); (5) Carilli & Barthel (1996); (6) 
Nagar et al. (2005); (6) De Ruiter et al. (1986); (7) Pellegrini (2005b); (8) Allen et al. (2006); (9) Allen et al. (2001); (10) Fabian et al. (2002); (11) Terashima 
et al. (2002); (12) Di Matteo et al. (2003); (13) Macchetto et al. (1997); (14) Bicknell and Begelman (1996); (15) Loewenstein et al. (2001); (16) Jones et al. 
(2007); (17) Taylor et al. (2006); (18) Giovannini et al. (1998); (19) Di Matteo et al. (2001); (20) Johnstone et al. (2005); (21) Zhao et al. (1993); (22) Blanton 
et al. (2003); (23) Zirbel & Baum (1995); (24); Simpson et al. (1996); (25) Evans et al. (2006) 



one measurement of the kinetic power was found, we have used 
their logarithmic average for our study, and increased the uncer- 
tainty accordingly. These are reported in Table|2] where a summary 
of the sample adopted is presented. For all these objects we define 
the Bondi power as 

Psondi = O.lAfBondiC^, (1) 

where the Bondi accretion rate is calculated from measures of gas 
temperatures and density under the assumption of spherical sym- 
metry and negligible angular momentum (see Allen et al. 2006, 
and references therein). For the nine AGN studied in Allen et al. 
(2006) the Bondi rate is calculated extrapolating the observed inner 
density profile of the X-ray emitting gas down to the Bondi radius. 
For the others, more luminous (and more distant) objects from the 
Rafferty et al. (2006) sample, the Chandra resolution corresponds 
to a size several orders of magnitude larger than the true Bondi ra- 
dius. For these objects, the Bondi power is estimated by assuming a 
inner density profile. This inevitably implies very large uncer- 
tainties on Peondi, that we have taken into account by increasing 
the nominal error bars associated with the measurements found in 
the literature. 

As mentioned in the introduction, our sample provides us with 
information on both the inner and the outer boundary condition for 
the accretion flow, as well as a reliable inventory of the energy out- 
put, both in radiative and kinetic form. In principle, if all these ob- 
jects were to be described by the same physical model for the ac- 
cretion flow (coupled to the jet/outflow), one should expect simple 
scaling relationships between the mass supply at the outer bound- 
ary (Bondi power), and the energy emitted in kinetic or radiative 
form from the accretion flow (Lkin and Lx, respectively). 



We first perform a partial correlation analysis in order to test 
which, if any, correlation between these quantities is statistically 
significant in our sample. We choose the partial Kendall's r cor- 
relation test for censored data sets dAkritas & Siebert 1996i The 
results of the correlation analysis are shown in Table |2] The cor- 
relations between the kinetic power and both Bondi power and 
X-ray luminosities are significant at more than 3-sigma level. In- 
terestingly, the correlation between Lx and Psondi does not ap- 
pear to be statistically significant once the common dependence of 
these two quantities on Lkin is accounted for. This is consistent 
with previous studies ( [Pellegrini 2005a[ ISoria et al. 2006) which 
failed to detect any clear correlation between Lx and Paondi, and 
might indicate that the X-ray luminosity, and thus the accretion 
power released radiatively, is not sensitive to the outer boundary 
condition, while it depends more strongly on the mechanical out- 
put of the flow. This is indeed expected, for example in adiabatic 
inflow-outflow (ADIOS) models ( |Blandford & Begelman T999) 
[Blandford & Begelman 2004) , but is also true for more general 
disk-wind models provided that a negligible fraction of the bind- 
ing energy of the accreting gas is converted directly into radia- 
tion dMerloni & Fabian 2002llKuncic & Bicknell 20"04l l, as we will 
show in more detail in section |4] 



On physical grounds, one might expect that both the black 
hole mass and the accretion rate should determine the power chan- 
neled through the jet, thus enforcing the need to perform mul- 
tivariate statistical analysis on any multi-wavelength database of 
AGN, a point already discussed in MHD03. However, we notice 
here that the sample in question spans only a very limited range 
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Table 2. Results of partial correlation analysis 





Variables 








Correlation 


X 


Y 


Z 




T 




Pnull 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


Log Lx 


Log Z/Kin 


Log 


PBondi 


0.49 


0.1818 


7.0 X 10-^ 


Log Psondi 


Log Lkiii 


Log 




0.5872 


0.2815 


3.7 X 10^2 


LogLx 


Log PBondi 


Log 




0.1289 


0.1422 


0.3647 



NOTE: Col. (1): Variable X. Col. (2): Variable Y. Col (3): Variable Z. Correlation between variables X and Y is studied, taking into account the mutual 
correlation of X and Y with Z. Col. (4)-(6): Results of partial correlation analysis, giving the partial Kendall's r con'elation coefficient, the square root of the 
calculated variance cr/f . and the associated probability P^uU for accepting the null hypothesis that there is no correlation between X and Y. 



of black hole masse|3- Attempts to quantitatively account for the 
mass dependence of the kinetic luminosity based on such a sam- 
ple were made, and yielded negative results. Therefore in what 
follows we limit ourselves to a simple bivariate correlation anal- 
ysis between Eddington-scaled quantities. Future, more detailed, 
studies of the physical connection between jet kinetic power and 
black hole accretion processes will greatly benefit from samples 
with much wider distributions of black hole masses. 

Figure [T] shows the relationship between Bondi power and ki- 
netic luminosity, both in units of the Eddington luminosity, I/Edd ~ 
AnMBHiTipc/aT, where rUp is the mass of a proton and ctt is the 
Thomson cross-section. The best fit linear regression slope, esti- 
mated via a symmetric algorithm that takes into account errors on 
both variables (see MHD03) gives: 

log(Lkin/iEdd) = (1.6i;^;^)l0g(PBondi/iEdd) + {1 -itU) (2) 

consistent, within the l-cr uncertainty, with the best fit slope found 
by Allen et al. (2006) (1.3lg;27). 

We have then studied the relationship between kinetic power 
and bolometric (radiated) luminosity. The bolometric correction for 
objects of such a low power is not well known. In general terms, it 
has been argued that low luminosity AGN in nearby galaxies dis- 
play a spectral energy distribution markedly different from clas- 
sical quasars, lacking clear signs of thermal UV emission usually 
associated to standard, optically thick accretion discs (see e.g. Ho 
2005, and references therein; for a contrasting view, Maoz 2007). 
The bolometric luminosity is probably dominated by the hard X- 
ray emission, but quantitative assessments of the bolometric cor- 
rections for low luminosity AGN based on well defined samples 
are still missing. Here, for the sake of simplicity, we adopt a com- 
mon 2-10 keV bolometric correction factor of 5 for all objects in 
our sample, and define Ax = SLx/^Edd- 

The best fit linear regression slope, estimated via a symmetric 
algorithm that takes into account errors on both variable|fl gives 
the following result: 

log(Lkin/LEdd) = (0.49 ± 0.07) log Ax - (0.78 ± 0.36) (3) 
with an intrinsic scatter of about 0.39 dex. 

^ This is likely the result of a specific selection bias, as the method used 
to estimate the average jet kinetic power can only be used effectively for 
radio galaxies at the center of bright X-ray emitting atmospheres, as those 
of clusters of galaxies (X-ray surface brightness selection). This in turn tend 
to select the more massive elliptical galaxies, with large central black holes. 
^ We assume here that the errors on Ax ai'e dominated by the statistical 
uncertainties on the black hole mass determinations, which we estimate to 
be of the order of ^^0.2 dex from the intrinsic scatter in the M-cr relation, 
see Tremaine et al. (2002). 




-5 -4 -3 -2 -1 

Log Peondi/LEdd 



Figure 1. Eddington-scaled kinetic power vs. Bondi power for the AGN in 
the sample. The thick solid line shows the best fit linear regression, with 
one-sigma uncertainties (thin solid lines). Red open circles with larger eiTor 
bars indicate those objects for which the Bondi radius is much smaller than 
the resolution and Bondi power is calculated by extrapolating inwards the 
observed density assuming a profile. 



The best fit relationship together with the data is shown in Fig- 
ure |2] and is inconsistent with a linear relationship at more than 
4-a level. Thus, not only are these SMBH more radio loud the 
lower their accretion rate (a fact already suggested by many previ- 
ous works both on nearby AGN and X-ray binaries: Ho & Ulvestad 
2001; Gallo, Fender & Pooley 2003; MHD03; Panessa et al. 2007), 
but their "kinetic loudness", i.e. the ratio of kinetic to bolometric 
powers (akin to a jet to accretion efficiency) is also a decreasing 
function of the dimensionless accretion rate onto the black hole. 
Even more interestingly, the best-fit slope coincides with the the- 
oretical expectations of radiatively inefficient, "jet-dominated" ac- 
cretion modes ( Fender, Gallo & Jonker 2003) . We will come back 
to this point in greater detail in section|4] 

As discussed in the introduction, many recent works have 
explored scaling relationships between black holes of stellar 
mass and supermassive ones by studying multidimensional cor- 
relation among different observables (MHD03; McHardy et al. 
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Log Ax 

Figure 2. Eddington-scaled kinetic power vs. bolometric nuclear luminosity 
for the AGN in the sample (black solid circles) and for the BHXRB Cyg X- 
1 (empty red star). The solid line shows the best fit linear regression, while 
the dashed line is the L^in = 5Lx relation shown as a term of comparison. 



2006), supporting the notion that AGN are indeed scaled-up 
galactic black holes. However, physical models for the disc- 
jet coupling in BHXRB based on the observed correlations be- 
tween radio and X-ray luminosity ( [Fender, Gallo & Jonker 2003] 
[Kording, Fender & Migliari 200 6) all face a large uncertainty due 
to the lack of reliable measurements of the jet kinetic power. In 
this context, it is interesting to include in Figure |2] the only GBH 
for which a measurement of the kinetic output has been made, Cyg 
X-1 dGallo et al. 20051 1, which turns out to be consistent with the 
relationship derived from the AGN sample. Clearly, systematic ef- 
forts to estimate kinetic power of BHXRB in the low/hard state are 
needed in order to assess their similarity with radio mode AGN. 



3 THE RELATIONSHIP BETWEEN KINETIC POWER 
AND RADIO CORE LUMINOSITY 

If the above correlation ([3) directly reveals fundamental physical 
properties of jet-producing AGN of low power, it still shows a non- 
negligible intrinsic scatter. On the other hand, one should expect 
a more direct relationship between the nuclear radio core emis- 
sion and the larger scale kinetic power, as both originate from 
the jet. All theoretical models for AGN flat-spectrum compact 
jet cores ([Blandford & Konigl 19791 IFalcke & Biermann 19961 
[Heinz & Sunyaev 2003) predict a dependence of the radio luminos- 

17/12 

ity on the jet power in the form Lr oc L^.^^ . The current sample 
provides by far the best opportunity to test these predictions. 

A Kendall's tau correlation test reveals that the kinetic power 
is correlated with the observed radio core luminosity LR.obs, with 



-Pnuii = 9.2 X 10~^ (see the empty circles in Figure We 
have fitted the data with a linear relationship: log Lkin = ^obs + 
-Bobs log //R.obs, once again making use of a symmetric regression 
algorithm that takes into account errors on both variables. We ob- 



tain Aobs = (22.1 ± 3.5), Bobs = (0.54 ± 0.09), with a large 
intrinsic scatter of cr = 0.47. 

Such a correlation, however, must be at some level biased by 
relativistic Doppler boosting of the radio emission in the relativistic 
jets. An alternative way to proceed would be to use indirect estima- 
tors of the nuclear radio core luminosity which are less affected by 
relativistic beaming dHeinz & Grimm 2005 1 , as, for example, the 
multivariate relation between BH mass, radio core and hard X-ray 
luminosity, the so-called 'fundamental plane' (FP) of active black 
holes (MHD03). Recent analysis of this correlation (Heinz & Mer- 
loni 2004; Kording, Falcke & Corbel 2006 [KFC06]; Merloni et al. 
2006) have shown that both Doppler boosting and sample selection 
play a crucial role in the exact determination of the intrinsic corre- 
lation coefficients of the FP, which also need to be accounted for. 
In the Appendix, we discuss in detail a possible way to overcome 
such a bias with the aid of a Monte Carlo simulation of the sam- 
ples used to derive the FP relation. That study allows us to estimate 
statistically the intrinsic (un-boosted) radio core luminosity of the 
AGN jets as a function of their (mean) Lorentz factor, Fm in a way 
that can be approximated by the following expression: 



log Lr,] 



(1-0.14 log Fm) [Crx log Lx + Crm log A^BH 

+CR(rm), 



(4) 



where Lr,pp is the intrinsic (un-boosted) radio core luminosity of 
the jet at 5 GHz, Lx the nuclear 2-10 keV intrinsic (un-absorbed) 
luminosity and Fm the mean Lorentz factor of the jets. 

In what follows, we will adopt the specific version of the FP 
relation derived from a sample of low luminosity AGN only, i.e. 
free from the bias introduced by the inclusion of bright, radiatively 
efficient AGN or QSOs (see discussion in KFC06). For that, the 
correlation coefficients are ^rx = 0.71, ^rm = 0.62, slightly 
different (but only at the 1-ct level) from those found in MHD03. 

Given Eq. l|4j, assuming a distribution of Lorentz factors for 
the AGN jets (or its mean, provided that the distribution is not too 
broad), we determine the "true" relationship between the Kinetic 
luminosity and radio core luminosity by fitting the 15 data points 
in our sample with the linear relationship 



log Lkin = Aint(rni) + -Bint(rm) log I/R,FP. 



(5) 



The fitted values for the intrinsic slope. Bint, as a function of Fm, 
are shown as a dot-dashed line in the bottom panel of Figure [3] 
From this we can see that the higher the mean Lorentz factor of 
the jets, the steeper must the intrinsic correlation between kinetic 
power and jet core luminosity be, and the larger the discrepancy 
with the measured slope, Bobs = 0.54 ± 0.09 (solid lines in Fig- 
ure O, obtained using simply the observed radio core luminosity, 
without any attempt to correct for relativistic beaming. In fact, such 
a discrepancy between the intrinsic and the observed slopes of the 
Lr - Lkin relation is indeed expected if the 15 sources of our sam- 
ple harbor relativistic jet randomly oriented with respect to the line 
of sigh0. 

In order to show this quantitatively, we have simulated (10'' 
times) the observed sample, assuming an underlying relationship 
between intrinsic radio core luminosity and kinetic luminosity 
given by Eq. l|5). In order to do that we have assumed distribu- 
tion distances, black hole masses radio and X-ray core flux limits 
that closely resemble the observed ones. We then Doppler-boost 



* It is worth noting here that the sample has no a priori selection against 
beamed objects (indeed, 3C 84, the radio source at the core of NGC 1275, 
is most likely beamed. Kirchbaum et al. 1992). 




Figure 3. Bottom panel: probability density of the observed slopes of the 
^Kin - ^R,obs relation (colored contours): For each value of the mean jets 
Lorentz factor, Fm, the intrinsic correlation has a slope -Bint determined 
from the best estimate of the beaming correction of the fundamental plane 
relation used as estimator of the intrinsic radio core luminosity (see text for 
details). The dot-dashed line in the lower panel shows this Bint, while the 
solid line shows the observed slope of the real data set with its uncertainty 
(dashed lines). The upper panels show then the integrated probability of 
observing 0.45 < B^ha < 0.63, given Bint andPni. 

the intrinsic radio core luminosity picking F from a normal dis- 
tribution with mean Fm and variance ar ~ O.lFm- Fits of the 
^kin-iR.obs Correlation for the 10* simulated samples result in a 
distribution of slopes Bobs as a function of Fm (shaded areas in 
the lower panels figure [3l(. We can now assess the probability of 
observing Bobs = (0.54 ± 0.09) for any value of Fm, by simply 
integrating these (properly normalized) distributions in the range 
0.54±0.09. The result of such an integration is shown in the upper 
panel of Figure |3] 

The distribution of the simulated sources in the L^^obs/ L-r. 
vs. Lkin plane shows tantalizing evidence of a discrepancy with the 
observed sample, with a slight deficit of luminous, de-boosted (i.e. 
seen at large angles with respect to the line of sight) objects. Sta- 
tistically, we found that this discrepancy mainly effects the normal- 
ization of the intrinsic Lj^in-I/R relation, rather than its slope. The 
reason for such a discrepancy is not clear at this point, but we be- 
lieve it may signal a problem in the modellization of the selection 
criteria of the sample, rather than a problem with the FP scaling. 
This is by no means surprising, in particular given the lack of any 
realistic constraint on the distribution and selection functions of the 
kinetic power measurements. 

From this Monte Carlo test, we can reach the following con- 
clusions: 

(i) For any given Fm, the beaming-corrected FP relation l|4) can 
be used to derive the intrinsic radio core luminosity of the jets, 
Lr,fp', 

(ii) The correlation between Li^in and I/r,fp in our sample has 
a slope i3int(rm), that is an increasing function of the mean jet 
Lorentz factor; 

(iii) The relation between Lkin and observed radio core luminosity 
is flattened with respect to the intrinsic one (i.e. Bobs < Bint) be- 
cause of the unaccounted-for Doppler boosting in the small sample 
at hand. We have corrected for such an effect statistically: Figure[3] 




Log Lr [erg s i] 

Figure 4. Kinetic power, as measured from estimates of pV work by the 
jet on its suiTounding hot, X-ray emitting, atmosphere versus radio core lu- 
minosity of the jets at 5 GHz. Red open circles show the observed radio 
luminosity (not-corrected for relativistic beaming), while the blue filled cir- 
cles refer to the radio luminosity estimated through the "fundamental plane" 
relation derived from the KFC sample, and corrected for relativistic beam- 
ing, with a mean Lorentz factor of the jets of Fm ~ 7 (see text for details). 
Thick (and thin) solid lines show the best fit to this latter relationship, with 
1-CT uncertainty; dashed lines show the best fit linear relation when the ob- 
served radio luminosity is considered. The value of the intrinsic scatter for 
the two cases is reported. 

shows that a fully consistent interpretation of the data can be given 
for a broad range of mean Lorentz factors, with a broad probability 
distribution that peaks at Fm — 7, corresponding to the following 
intrinsic correlation: 

log Lkin = (0.81 ±0.11) log Lr + ll.gtti (6) 

Reassuringly, the probability distribution of mean Lorentz factors, 
although broad, clearly excludes low values of Fm- More detailed 
studies of the relativistic speeds of Blazars emitting regions give 
similar results (see e.g. Cohen et al. 2007). 

Interestingly, the relationship between kinetic power and ra- 
dio luminosity gets tighter if one uses the FP relation as an estima- 
tor of the un-boosted radio flux, as shown in Figure |4l the intrin- 
sic scatter is reduced from 0.47 dex to 0.37 dex. This is expected 
if this steeper correlation reflects more directly the intrinsic rela- 
tion between jet kinetic and radiative power, and lends support to 
our statistical method to correct for the bias induced by relativistic 
beaming. 

Finally, we note that the slope of the correlation between 
Lr.pp and Lj^in is consistent with the predictions of syn- 
chrotron models for the flat spectrum jet cores according to which 
( |Heinz & Sunyaev 2003| l Lr oc /^^i7+8«^)/i2j^,/-a,^ ^j^^^.^ 
is the radio spectral index. For flat spectrum cores, like the ones 
we are considering here, ~ and the logarithmic slope of the 
Lkin ~ Lr correlation is expected to be 12/17 — 0.71, within 
1-(T from our estimate. 
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The intrinsic relation between radio and kinetic luminosity of 
the jets can also be used to determine the total radiative efficiency 
of the jet cores, 77jct, the ratio of radiated to kinetic luminosity. This 
will of course depend on the maximal frequency, i^max, up to which 
the flat radio core spectrum extends. Assuming a constant spectral 
index (i.e. neglecting spectral steepening ati/K i^max) we have: 

logr?jot = (a, + 1) log(i.n,ax/5GHz) + 0.19 log Lr - 11.9 (7) 

Many nearby low-luminosity radio galaxies observed with the HST 
show unresolved optical nuclei the flux of which clearly correlates 
with the radio one ( Chiab erge, Capetti & Celotti 1999| l, suggesting 
that in those object the jet emission extends up to the optical band, 
and the corresponding radiative efficiency can be as high as 10- 
20%. This result is not substantially modified if we allow at least a 
non negligible fraction of the X-ray core luminosity to originate in 
the jet, as recently argued on the basis of either X-ray spectroscopy 
alone JEvans et al. 2006 ) or on model fitting of multiwavelength 
SED jWu, Yuan & Cao 2007[ >. This is because the available con- 
straints on broad band SED of low-luminosity radio AGN (mostly 
FR I) based on radio-to-optical-to-X-ray indexes (see e.g. Hardcas- 
tle & Worrall 2000; Heinz 2004; Balmaverde, Capetti & Grandi 
2006) show that the spectra must be concave everywhere between 
the radio and X-ray bands, and this alone ensures that the radiative 
output of the jet cores is dominated by the emission at the peak 
frequency, t'max, thus justifying the use of the above equation. 



4 DISCUSSION 

4.1 Testing accretion flow models 

In section |2] we have shown that there is a significant correlation 
between kinetic power and Bondi power, as well as between X-ray 
(bolometric) luminosity and kinetic power, while the Bondi power 
and the X-ray luminosity are not correlated. Here we try to interpret 
this fact, as well as the observed slopes of the statistically signifi- 
cant correlations within simple models for sub-Eddington accretion 
flows. 

The idea that low-luminosity AGN accrete in radiatively in- 
efficient fashion was first put forward by Rees et al. (1982) and 
Fabian & Rees (1995), while efforts were made at the same time 
to accommodate advection dominated accretion flows (ADAF; 
Narayan & Yi 1995) within the general framework of accretion 
theory. Very early on it was also recognized that such kind of adi- 
abatic (i.e. non radiative) accretion flows are prone to the pro- 
duction of strong outflows, so that in fact only a tiny fraction 
of the gas entering the black hole sphere of influence finally 
makes its way onto the black hole (Bla ndford & Begelman 1999] 
[Blandford & Begelman 20041. While the observational evidence 
for low radiative efficiency in many nearby galactic nuclei has 
grown significantly in recent years (see e.g. Ho 2002; MHD03; 
Pellegrini 2005; Chiaberge et al. 2005) there remains substantial 
uncertainty on the true dynamical state of such flows. 

Here we want to exploit the unique possibility offered by this 
sample with its comprehensive information on both inner and outer 
boundary conditions for the accretion flow, as well as its reliable in- 
ventory of the energy output. We start from a very general accretion 
disc model coupled to a disk wind such as those studied in detail 
in Kuncic & Bicknell (2004). There the interested reader will find 
an analytic description of the most general MHD disk accretion, 
specifically addressing the relationship between radial and vertical 
mean field transport of mass, momentum and energy. 



In what follows, we will use as a starting point Eq. (84) of 
Kuncic & Bicknell (2004) for the rate of energy released by the 
viscous torques within the disc as: 



3GMbhMou 



<K{r) 



3»K 



M^ir')vKir') 



dr'(8) 



where Mout is the accretion rate at the outer boundary, = 1 — 
{rin/ry^'^, Ok is the Keplerian angular velocity D,k = vk/t = 
^GMIr\ Mn,{r) is the mass loss rate in the outflow, and we 
have neglected any additional torque at the disc inner boundary, as 
well as all additional terms due to vertical Poynting flux and re- 
irradiation of the disc from outside, for the sake of simplicity. 

The first term on the right hand side of Eq. ^ represents the 
binding energy flux associated with the net mass flux, and is famil- 
iar from standard accretion disc theory ( [Shakura & Sunyaev 1973[ l, 
while the second term is related to the flux of kinetic and gravita- 
tional energy in the outflow. 

We adopt here the simplest phenomenological model for the 
accretion rate through the disc: 



Ma{r) 



Mout, 
Mout ( 



r > Tout 
r < Tout 



(9) 



with < p < 1 piandford & Begelman 1999^ , where rout rep- 
resents a critical radius beyond which mass outflow is negligible. 
Correspondingly, mass conservation implies that the mass outflow 
as a function of radius can be expressed as = Afout — Ma = 
Afout [1 — {r /rout ■ With such a simple ansatz, the kinetic energy 
properties of the outflow closely resemble those of the adiabatic 
ADIOS solutions. 

By integrating Eq. ([Sj from the inner radius, rin to rout, we 
obtain for the power dissipated within the disc: 



rout 



GMMou 



dr 
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(l-p) (2p + l)(l-p) 
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(10) 



where we have introduced the parameter ^ = In the above 

expression. Pa represents the total power released by the accretion 
process, while Lmcch is the mechanical energy carried by the out- 
flow. The latter is given by 

— / -K(r)r— dr 
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(11) 



where Voo is the ratio of the terminal wind speed to the Keplerian 
velocity, here assumed to be independent of r. 

The radiative power depends on micro-physical properties of 
the disc (capability of heating electrons as opposed ions), on its ver- 
tical structure and on the radiative processes involved. Radiatively 
efficient discs, as those discussed in Kuncic & Bicknell (2004), 
would obviously have Lboi = Pd, and we can already guess that 
this option is not allowed by the observed correlations. Indeed, in 
the limit of very large rout, C ^ 1, we have that both the me- 
chanical luminosity of the disc, Lmcch and Pd scale as the accre- 
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tion rate through the disc in its inner part Afa ~ Mout^*', so that 
rather than the observed Lmcch oc L^'^i (see sec- 
tion|2]l. The only way in which this can be altered is by considering 
the effect of magnetic fields in the accretion flow. As it was pro- 
posed by Merloni & Fabian (2002) and discussed also in Kuncic & 
Bicknell (2007), strong magnetic field threading a thin disc may re- 
move energy and angular momentum (and very little mass) without 
radiating efficiently, and in so doing generating powerful, Poynting 
dominated outflows. However, the specific scaling between kinetic 
and radiative power in this case depends in a complicated way on 
the dynamo properties of the disc, and is not yet as fully understood 
as the adiabatic (ADIOS) solutions we discuss here. 

This suggests a simple phenomenological scaling for the ra- 
diative output of mass-losing low luminosity discs, that can be writ- 
ten as log (Lboi/iEdd) = C+7log(Pd/iEdd)- The normalization 
C is related to the value of the dimensionless accretion rate above 
which radiatively inefficient accretion ceases to be a viable solution 
due to the increase of disc density and cooling rate. The value of 
the parameter 7, instead, depends on the radiative cooling proper- 
ties of the accreting gas, and within our approach has to be fixed by 
comparison with observations. 

To fully characterize the model, and make it suitable to be 
tested against the available data, one still needs to specify the rela- 
tionship between rout and Alout ~ AfBondi. In adiabatic solutions 
with outflows, the radius beyond which radiative losses are signif- 
icant and the adiabatic self-similar solution breaks down will de- 
pend on the outer supply of mass and angular momentum and on the 
cooling properties of the outer disc ( [Blandford & Begelman 2004] l. 
In our case, we simply assume: 



log^ = A + Blog(Mout/A/Edd) 



(12) 



and use the data to put limits on A and B. 

Our accretion/outflow model has the desired properties of sim- 
plicity and generality, which make it suitable to be tested against the 
dataset at hand. Indeed, variations of the indexes p and 7 allow us to 
span the full range of dynamical accretion models from pure ADAF 
(p = and 7 ~ 2) to ADIOS (1 > p > and 7 ~ 2) to radiatively 
efficient, wind dominated discs(l>p>0 and 7 ~ 1) down to 
standard, radiatively efficient, conservative thin discs (p — and 
7 = 1). The parameters A and B, in turn, give us an handle on 
the large-scale geometrical properties of the flows, with C repre- 
senting an overall normalization of the radiative efficiency of the 
system. The sixth and final parameter of the model, Voc, turns out 
to be strongly degenerate with the normalization parameters A and 
C. Given the amount of information in the data, we have chosen to 
keep its value fixed, either to a value of 1 (terminal velocity equal 
to the local Keplerian speed) or 2 (terminal velocity twice the local 
Keplerian speed), and perform two separate fits. Detailed dynami- 
cal modelling of the jets and outflows at different luminosities and 
black hole masses will be needed in order to constrain its value 
reliably. 

In summary, our model has just five free parameters: p, 7, 
C, A and B. We have fitted this model to the observed relation- 
ships between kinetic power and Bondi power, and between kinetic 
power and bolometric luminosity. In Figure |5] we show the best fit 
to the data for the two independent correlations with the model de- 
scribed above (in the case Voc ~ 2, blue solid line and Voo ~ 2, 
green solid line). Our simple illustrative model provides an excel- 
lent description of the data both in the case of Voo = 1 and Voc ~ 2. 
The best-fit values for the model parameters (with their 90% confi- 
dence errors, calculated by letting each single parameter vary while 
fixing all others to their best fit value) are given by: p — 0.85^q :tg. 
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Figure 5. The best fit accretion/outflow models (blue and green solid line 
for tioo = 2 and 1, respectively) are plotted alongside the simple best fit 
linear relations among Eddington-scaled kinetic power, Bondi power and 
bolometric luminosity (see text for details). 



A = 0.13 ± 0.30 B = 0.63 ± 0.10 and 7 = 1.98 ± 0.06 and 
C = 2.27 ± 0.24 in the case = 2, and p = 0.92l° :[^, 
A = 0.1 ± 0.19 B = 0.45 ± 0.10 and 7 = 2.04 ± 0.06 and 
C = 2.00 ± 0.22 in the case Va^ = 1. 

A crucial prediction of the model is the dependency of rout on 
the outer accretion rate. According to our fit, we have, in units of 
the Schwarzschild radius rs = 2GM/(?, 

-0.63 



rout 



6 X 10" 



"tout 

10-4 



(13) 



where rhout = Afout/MEdd. For the objects in our sample, this 
is reassuringly smaller than either the inferred Bondi radius and 
the radius beyond which radiative losses are significant and adi- 
abatic self-similar outflow dominated solutions are not viable. In 
the sample under study, the two brightest objects (Cyg A and Hy- 
dra A) have rhout ~ 10^^'^ and therefore outer outflow radii of 
just a few tens of Schwarzschild radii. They may therefore qualify 
as "intermediate" objects in a transition between radiatively effi- 
cient (hidden QSO-like) and inefficient AGN. In this case, the self- 
similar assumption becomes poorer, and clear signs of the pres- 
ence of an outer standard disc may be expected ( [Ogle et al. 1997[ 
ISambruna et al. 2000t . 

As we have discussed above, a clear distinction exists between 
those parameters that describe the nature of the outer boundary con- 
ditions of the flow and solely determine slope and normalization of 
the I/Kin-fBondi relation {A, B and p) and those that instead de- 
scribe the nature and radiative properties of the inner accretion flow 
and determine the slope of the Lx - I/Kin relation (7 and p). This 
might explain why we found no statistically significant correlation 
between radiative output and Bondi power. 

It is worth stressing here that the above model represents a 
very schematic description of the physical properties of a low- 
luminosity accretion flow and the accompanying powerful outflow. 
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As mentioned above, a full exploration of all possible theoretical 
alternatives for such flows should also take into account dynami- 
cally important magnetic fields and their role in the jet production 
mechanisms. The reason we have chosen not to do so here is merely 
the quality of the data available, which would not have allowed to 
put meaningful constraints on the parameters of such complicated 
MHD disc-outflow solutions. What we have instead shown is that 
future, better samples like that assembled here promise to be ex- 
tremely useful for constraining the complex physical properties of 
accretion at the lowest rates. 



4.2 The kinetic power output of Sgr A* 

The quiescent radio and X-ray luminosities of the ~ 4 x 10'° Mq 
nuclear black hole at the center of the Milky Way, Sgr A* 
(logl/R ~ 32.5 and logLx — 33.3) would be both consistent 
with a Bondi accretion rate of about 10~® Mq /yr, in line with the 
Chandra measurements of the inner hot gas properties of the galac- 
tic center (Baganoff et al. 2001 ; Baganoff et al. 2003). If the same 
model outlined in section |4. 1 [ applied to Sgr A*, it would predict an 
outer radius of the outflow of about Tout ~ 2.5x10"^ Schwarzschild 
radii (~ 0.025" at the distance of the galactic center) and an in- 
ner accretion rate Ma(j-in) = A/outC'' 4 x 10"^ Mo/yr, 
also consistent with the gas density in the inner accretion flow 
implied by linear polarization measurements ( lAitken et al. 20001 
IBower et al. 2003l|Bower et al. 2005l >. Altogether, this would indi- 
cate that the SMBH at the galactic center is the source of about 
5 X 10^* ergs s^^ of mechanical power (similar to what estimated 
within jet model fits to Sgr A* SED, see e.g. Falcke & Biermann 
1999, Falcke & Markoff 2000), or about 1.5 supemovae every 10^ 
years. Such a mechanical power input into the galactic center could 
play a significant role in the production of the TeV 7-rays recently 
observed by the HESS (High Energy Stereoscopic System) collabo- 
ration dAharonian et al. 20041 [Atoyan & Dermer 2004[ >, as well as 
in the heating of the hot (> 8 keV) diffused plasma detected by 
Chandra JMuno et al. 20041 ). 



4.3 The effect of variability on the observed correlations 

In section [3] we have concentrated our attention on the assess- 
ment of the effects of relativistic beaming on the measured slope 
of the Lkin-^R correlation. However, one should expect a second 
source of scatter in any correlation between kinetic power and nu- 
clear luminosity (in any waveband). The large scale power is an 
average over the typical age of the cavities and bubbles observed 
in the X-ray atmosphere of the galaxies in our sample. Such an 
age may be assumed to be of the order of either the buoyant rise 
time, or the sound crossing time, or the refill time of the radio 
lobes dBirzan et al. 20041 . These estimates typically differ by about 
a factor of 2, and, for the sample in question, lie in the range 
^age ~ 10^ — 10* years. On the other hand, the core power is 
variable on time scale much shorter than that. The bias introduced 
by this fact can be quantified as follows. 

First of all, we notice that AGN X-ray lightcurves, where most 
of the accretion power emerges, show a characteristic rms-flux re- 
lation which implies they have a formally non-linear, exponential 
form ( [Uttley, McHardy & Vaugh an 2005 ), and the luminosity fol- 
lows a log-normal distribution (Nipoti & Binney 2005| ). Let us, for 
the sake of simplicity, define as L{t) the instantaneous luminosity 
of the AGN, be it bolometric, kinetic or radio, under the implicit as- 
sumption that very close to the central engine jet and accretion are 



strongly coupled, so that all of them follow a log-normal distribu- 
tion. Thus, log L is normally distributed with mean m and variance 
af . The measured kinetic power Lkin discussed in this work being a 
long-term time average, it should be determined by the mean of the 
log-normally distributed L, i.e. (L) — fi — exp (/if + erf /2). Be- 
cause the log-normal distribution is positively skewed, the mode m 
of the distribution of L, i.e. the most likely value of a measurement 
of it. Lobs, is smaller than the mean: Lobs — m — exp {fii — erf). 
Therefore the most likely value of the ratio of the mean to the ob- 
served luminosity is given by (Nipoti & Binney 2005; Uttley et al. 
2005): 

^-ei^' ^{aL. + lf', (14) 

where we have introduced the rms (fractional) variability of the ob- 
served lightcurve crfms — exp erf — 1, an easily measurable quan- 
tity. The higher the rms variability of a lightcurve, the higher is the 
probability that an instantaneous measurement of the luminosity 
yields a value smaller than (L) ( [Nipoti & Binney 2005) , and also 
the higher the most likely ratio between the two values. 

Before discussing what are the appropriate values of <j^^^ to 
be used in Eq. il4\ . we also note that the observed slope of the 
(L)-Lobs correlation may deviate from unity for a log-normal dis- 
tribution, thus skewing any observed correlation between mean and 
instantaneous power, like those we have discussed so far. It is easy 
to show that from Eq. l ll4t we obtain: 

aiog(L) ^ I 3 ""Ls ^logO-fms ^jj^ 
9 log Lobs 2(1+ CT?ms ) 9 log Lobs ' 

Log-normal AGN variability may thus skew the observed relation 
between jet average kinetic power and instantaneous core lumi- 
nosity much in the same way relativistic beaming does. Inspec- 
tion of Eqs. l ll4t and \151 suggests that the measured slopes of 
any correlation of the (average) kinetic vs, core power will be 
substantially affected by variability if, and only if, both crfms and 
d log crfjns/9 log Lobs are at least of the order of unity. 

Unfortunately, very little is known observationally about the 
variability amplitude of AGN on very long timescales, especially 
for AGN of low luminosity as those considered here. Brighter AGN 
(Seyferts), on shorter timescales (1-10 years), have indeed rms 
variability amplitudes that rise steeply with decreasing luminosity, 
from about afms ~ lO"'^ for Lx « 10** up to crfms ~ 10~^ 
for Lx ^ 10*^ dNandra et al. 19971 IMaAowitz & Edelson 200II ). 
However, no evidence is yet found of such a trend continuing 
down to lower luminosities. On the contrary, suggestions have been 
made that crf^jg may flattens out at lower Lx ) [Papadakis 2004[ 
IPaolillo et al. 20041 ) at values of a few times 10~ for the typical 
X-ray luminosity of the objects in our sample. This would imply 
that the observed correlations between large scale and core powers 
are not skewed by variability bias. 

However, we should also consider the possibility that, 
if measured on the tagc timescale, rather than on years or 
decades, AGN variability amplitude can be much higher, per- 
haps exceeding unity (Nipoti & Binney 2005). Specifically for 
the case of radio galaxies, various indirect pieces of evidence 
of intermittency have been presented, such as the ripples and 
shock waves detected in the X-ray emitting atmosphere of M87 
dForman et al. 2005) , or the number vs. size counts of small radio 
galaxies ( [Reynolds & Begelman 1997) , both obviously relevant for 
our current discussion. 

All AGN Power Spectral Densities (PSD: the variability 
power P{i^) at frequency u, or timescale l/v) are best fitted by 
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a power-law of index —1 (P{y) oc i^^^) on long timescales, 
which breaks to a steeper slope on timescales shorter than a break 
timescale, which itself appear to be correlated with the central black 
hole mass and accretion rate (tbrcak oc A/|jj/Lboi, McHardy et 
al. 2006). Indeed, if the PSD of AGN at very low frequencies, 
V ~ 1/tage, are a seamless extension of the flicker noise observed 
on days-years timescales, then their rms variability amplitude could 
be very large: (Trms ~ 1 — 10. 

If we define the duty cycle of L{t) as 5 = {L)'^ / {Li^) = 
(1 + cTrnis) ^ dClottl & Ostrlkcr 2001) , then large rms amplitudes 
would imply very bursty lightcurves, with very short duty cycles 
5 ~ 10~^ — 10~^. Within this picture, the measured Lkin would 
be dominated by very short, "quasar-like phases", rather than by 
quasi-continuous radio mode AGN activity. 

Although this possibility cannot be easily ruled out, we feel it 
is currently disfavored for two main reasons: first of all, it would 
imply a typical QSO lifetime for these objects of < 10^ yrs, 
much less than current estimates JMartini 2004b : secondly, if the 
variability properties of low luminosity AGN are scaled up ver- 
sions of those of low/hard state GBH, as it is the case for for bright 
(high/soft state) objects f McHardy et al. 2006) , then we should ex- 
pect a second, lower frequency break to white noise slope {P{v) = 
const.) in the PSD ( [Utfley, McHardy & Vaughan 2005} . Such a 
break should be present at timescales fflat between 10-100 times 
longer than tbrcak, but still much smaller that tagc- In this case, 
~ iV[2 + ln(tflat/tbrcak)] 6iV, whcre TV is the fre- 
quency independent amplitude of f x P(i') in the flicker noise 
part of the PSD, which has been measured for many AGN and 
GBH to be of the order of a few times 10~^ ( [Papadakis 2004[ 
IDone & Gierlihski 2005 1 . Thus, at most, we should expect (Trms ~ 
a few times lO"'^ also for low luminosity AGN on very long 
timescales, implying a much higher duty cycle (not much smaller 
than unity). 

To summarize, the effects of the AGN intermittency on the ob- 
served relations between time-averaged kinetic power and instan- 
taneous core luminosities will be comparable to that of beaming 
if AGN have very bursty lightcurve (duty cycles less than 10%) 
on timescales comparable to age of the bubbles over which the ki- 
netic power is estimated, and the overall burstiness of the lightcurve 
(its measured rms variability) increases with decreasing luminosity. 
Although this possibility cannot be ruled out, there are reasons to 
believe that the rms variability of the low luminosity AGN in our 
sample is smaller than unity. If this is the case, the difference be- 
tween instantaneous and average power for these objects is of the 
order of (or smaller than) the systematic uncertainties on the power 
itself, and the observed correlations between Lkin and nuclear lu- 
minosities should not be strongly affected by variability. More work 
is needed to investigate this issue, which is beyond the scope of this 
paper. 



5 CONCLUSIONS 

We have presented a statistical analysis of a sample of 15 nearby 
AGN for which the average kinetic power has been estimated from 
the study of the cavities and bubbles produced by the jets in the 
IGM. In particular, rather than focusing on the relationship between 
the kinetic power and the IGM physical state, as was done before, 
we have tried to derive the relationship between jet kinetic power 
and nuclear properties of the AGN, specifically in terms of their 
black hole masses, 2-10 keV and 5 GHz radio core luminosities. 
Following our analysis, we reach the following conclusions: 



• A clear relationship exists between Eddington-scaled kinetic 
power and bolometric luminosity, given by: log(Lkin/i'Edd) = 
(0.49±0.07) log(Lboi/iEdd)-(0.78±0.36). The measured slope 
suggests that these objects are in a radiatively inefficient accretion 
mode. 

• We confirm previous claims of a correlation between Bondi 
power (i.e. accretion rate at the Bondi radius) and Kinetic lumi- 
nosity of the jets. Interestingly, there is no statistically significant 
correlation between Bondi power and bolometric luminosity apart 
from that induced by their common dependence on Lkin- 

• The observed correlations are in very good agreement with 
theoretical predictions of adiabatic accretion models with strong 
outflows. We have also shown that meaningful constraints on some 
specific physical properties of such models can be placed by fitting 
the observed data set. 

• The available measures of the average jet power provide a 
very useful tool to assess, in a statistical sense, both the jet radia- 
tive efficiency and the effects of relativistic beaming on the ob- 
served AGN jet population. Combining information on the kinetic 
jet power with estimators of the un-beamed radio flux of a jet core 
(as, for example, the so-called fundamental plane of active black 
holes, MHD03), we are able to determine both the probability dis- 
tribution of the mean jets Lorentz factor, that peaks at F ~ 7, and 
the intrinsic relationship between kinetic and radio core luminos- 
ity, logLkin = (0.81 ± 0.11) logLR + 11.9l:4;4 which is in good 
agreement with theoretical predictions of synchrotron jet models: 
ivkin C>C Where is the radio spectral index of the 
jet core. 

• The total radiative efficiency of the jet can be expressed as a 
function of the observed 5 GHz luminosity (Lr) and spectral index, 
a, as: r)i.t ^ 3.2 x 10-5(LR/1033erg/s)°-i^(i/^ax/5GHz)i+"-, 
where i^max is the high turnover frequency of the synchrotron emis- 
sion 

• The relation between Lkin and observed radio core luminos- 
ity is flattened with respect to the intrinsic one because of the 
unaccounted-for Doppler boosting in the small sample at hand. Re- 
sults of previous works that have studied the relationship between 
total extended radio power and kinetic power dBirzan et al. 20041 
IBirzan et al. 2006l l. found both flat slopes and very large intrinsic 
scatter (much larger than what we found here). As Birzan et al. 
(2006) noticed, such a scatter is "intrinsic to the radio data, for rea- 
sons that include radio aging and adiabatic expansion." This might 
be considered as a further argument is support of our approach of 
using core radio luminosities, instead, as estimator of the AGN ki- 
netic power. 

• The intrinsic variability of the AGN could affect our results 
due to the fact that the measured kinetic power is a long-term aver- 
age of the instantaneous power, and that the luminosity of accreting 
black holes are log-normally distributed, and therefore positively 
skewed. Very little is known about the long-term (10^-10** years) 
variability properties of low luminosity AGN. However, available 
estimates of rms variability as a function of luminosity on shorter 
timescales, scaling relations with galactic black holes as well as 
current estimates of quasar lifetimes all seem to suggest that the 
variability cannot be so extreme to affect the results of a correla- 
tion analysis more dramatically than relativistic beaming does. 

With the aid of these findings, quantitative assessments of kinetic 
feedback from supermassive black holes in the radio mode (i.e. at 
low dimensionless accretion rates) will be possible based on accu- 
rate determinations of the central engine properties alone, such as 
mass, radio core and/or X-ray luminosity. This will provide useful 
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constraints for AGN feedback models in a cosmological context 
(Merloni and Heinz 2007, in prep.). 
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APPENDIX A: THE EFFECT OF RELATIVISTIC 
BEAMING ON THE FUNDAMENTAL PLANE RELATION 
FOR ACTIVE BLACK HOLES 

We are concerned here with the following question: What is the 
effect of relativistic beaming on the determination of the slope of 
the fundamental plane (FP, MHD03) relation? Or, put it in another 
way, can we safely assume that the FP relation (in any of its in- 



carnations) is free from biases due to relativistic beaming, and thus 
suitable to use as a calibrator to estimate the intrinsic radio core 
luminosity? 

As originally discussed in MHD03, we introduce the "ob- 
served" FP relation as 

log I/R.obs = Crx log Lx + Crm log Mbh + ba, (Al) 

where the observed correlation coefficients depend, at the 1-sigma 
level, on the specific choice of sample, as discussed in KFC06. 

The strongest effect on the slope of any intrinsic correlation 
(i.e. any correlation between the intrinsic, un-boosted luminosities 
of jetted sources) comes from the inclusion/exclusion of object with 
jet axis close to our line of sight (see e.g. Heinz and Merloni 2004). 
If we define a "cut-off" angle, 9c, such that all object with inclina- 
tion 9 < 9c are excluded from a sample, then the slope of the ob- 
served correlation deviates from the intrinsic one as a strong func- 
tion of 6c once the Lorentz factor F > 1 /9c. 

In order to assess this effect, we have simulated a sample of 
nearby SMBH according to the selection criteria of MHD03 and 
Kording, Falcke and Corbel (2006) [KFC06], with a simple Monte 
Carlo routine. In our fiducial calculation, an effective cut-off an- 
gle of 9c — 15° is assumed, as these samples have been "cleaned" 
of relativistically boosted radio cores by excluding (or correcting 
for) all known Bl Lac objects. The radio fluxes are drawn from 
the observed local (z=0) flat spectrum radio luminosity function 
( [Dunlop & Peaco ck 1990). Distances are drawn from a distribution 
that closely resembles that of the soirees in our sample (see Ta- 
ble |2](, and objects with a flux smaller than 0.1 mjy at 5 GHz and 
10"^* ergs s~^ cm~^ in the 2-10 keV band are excluded. The ra- 
dio luminosity is then un-boosted assuming random orientation of 
the jet between 9c and 90°. The jet Lorentz factor is assumed to be 
normally distributed around Fm with variance or = O.lFm, and 
we studied the effects of varying Fm between 1.5 and 50. We fitted 
the simulated data with the (symmetric) OLS Bisector regression 
algorithm dlsobe et aL 19901 in order to determine the relationship 
between the intrinsic radio luminosity Lr and the observed one: 
log Lr = (fRR(rm) log Z/R,obs + feRR. As expected, we found that 
the larger the mean Lorenz factor of the objects, the more the ob- 
served correlations are skewed away from the ^rr = 1 intrinsic 
slope. In particular, the best fit slopes can be approximated by the 
following log-linear relation: 

Crr~ l-0.141ogFn„ (A2) 

that can be regarded as a "calibrator" of the fundamental plane re- 
lationship. 

The above relation l lA2t is plotted in both upper and lower 
panels of Fig. lAll as a solid line, while the dark red shaded ar- 
eas represent the 1-sigma contours. Figure IaTI also illustrates the 
effects of slightly varying our fiducial assumptions on the distri- 
bution of the Lorentz factors and on the cut-off angle. In partic- 
ular, the upper panel shows the results of increasing the width of 
the Lorentz factor distribution by a factor of 3 (from O.lFm to 
0.3Fm, dot-dashed lines). In this case the best fit slope is approx- 
imated by another log-linear with slope -0.16 (thick dot-dashed 
line). This translates into a difference in the FP slope from our 
fiducial case around Fm — 7 of just about 3%, much less than 
the statistical uncertainties. Similar is the case when different cut- 
off angles are considered. In the lower panel of Fig. lAll we show 
how the observed correlations are skewed due to relativistic beam- 
ing when 9c = 10° (dot-dashed contours and thick dot-dashed line) 
or 9c = 20° (dashed contours and thick dashed line). Also in these 
cases the difference with respect to our fiducial case corresponds to 
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Figure Al. The effects of relativistic beaming on the observed radio lumi- 
nosities in simulated samples resembling those used to derive the "funda- 
mental plane" relation. In the upper panel we show, as a function of the 
mean jet Lorentz factor, the best-fit slope "f the relation between in- 
trinsic and observed radio core luminosity for two values of the width of 
the r distribution: (Tp/F = 0.1 (solid contours and thick solid line) and 
or/r = 0.3 (dot-dashed contours and thick dot-dashed line). In the lower 
panel we show (,y{M. ^s. Fm for three values of the cut-off angle: dc = 10° 
(dot-dashed contours and thick dot-dashed line), dc = 15° (solid contours 
and thick solid line) and Be = 20° (dashed contours and thick dashed line). 

a difference in the "corrected" FP coefficients of less than ~3% at 

Vm ^ 7. 

In summary, for any value of mean jet Lorentz factor of the 
sampled AGN, Fm, we have shown that it is possible to derive sta- 
tistically an intrinsic (un-boosted) radio core luminosity based on 
the observed hard X-ray one, Lx and on the black hole mass, A/bh 
according to: 

logZ/R,FP = (1 — 0.14 log rm)KRX log Lx + Crm log Mbh] 

+CR(rm), 

where we have defined the new constant cr = 5rr,6r, + 6rr,. The 
term in the first parenthesis in the right hand side thus represents our 
simple way to estimate the relativistic beaming bias introduced in 
the samples used originally to define the FP relation. Its numerical 
value do indeed depends on the assumptions of our fiducial Monte 
Carlo model, but to an extent that is negligible when compared to 
the statistical uncertainties on the FP parameters themselves. This 
above relation, then allows a meaningful statistical test of the in- 
trinsic correlation between radio core luminosity and kinetic power 
of AGN jets, as we show in section[3] 



